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Abstract 

A Solar Sail Radiation Environment (SSRE) model has been developed for defining charged particle environments 
over an energy range from 0.01 keV to 1 MeV for hydrogen ions, helium ions, and electrons. The SSRE model 
provides the free field charged particle environment required for characterizing energy deposition per unit mass, 
charge deposition, and dose rate dependent conductivity processes required to evaluate radiation dose and internal 
(bulk) charging processes in the solar sail membrane in interplanetary space. Solar wind and energetic particle 
measurements from instruments aboard the Ulysses spacecraft in a solar, near-polar orbit provide the particle data 
over a range of heliospheric latitudes used to derive the environment that can be used for radiation and charging 
environments for both high inclination 0.5 AU Solar Polar Tmagpr mission and the 10 AU LI solar sail missions. 
This paper describes the techniques used to model comprehensive electron, proton, and helium spectra over the 
range of particle energies of significance to energy and charge deposition in thin (<25 /xm) solar sail materials. 

1*0 Introduction 


Materials exposed to space radiation environments degrade due to the accumulation of displacement damage and 
dissociation of chemical bonds produced by passage of charged particles through the material. Significant localized 
damage can also occur due to electrostatic discharging, arcing, and treeing processes when electric fields generated 
by accumulated charge densities exceed the dielectric strength of the material. Evaluation of radation susceptibility 
and risks of damage due to arcing for candidate solar sail materials require radiation environment models for 
establishing laboratory tests of radiation induced changes in material properties, modeling of radiation dose as a 
function of depth in materials, and modeling of the charging process. 

This paper describes development of the current version (Revision A) of the Solar Sail Radiation Environment 
(SSRE) model developed to characterize solar wind and energetic particle (electron, proton, and helium) 
contributions to mission accumulated dose and internal charging over an energy range O.OlkeV < E < 1 MeV. 
SSRE provides environments for Solar Polar Imager’s planned 0.5 astronomical unit (AU) solar polar orbit and the 
proposed ~1.0 AU first Lagrange point (LI) solar sail missions (LI Diamond and Solar Sentinel missions). 

2.0 Data and Model Construction 


Solar wind and energetic particle instruments aboard the Ulysses spacecraft provide the approximately 12.6 
years of solar wind and energetic charged particle data used to derive the radiation environment. Ulysses is the only 
spacecraft to sample high latitude solar wind environments far from the ecliptic plane and is therefore uniquely 
capable of providing the information necessary for defining radiation environments for the Solar Polar Imago- 
spacecraft. The Solar Wind Observations Over the Poles of the Sun (SWOOPS) instrument [Bame et al ., 1992] 
provides the low energy electron and ion component of the solar wind and the Low-Energy Magnetic Spectrometer 
(LEMS) and Low Energy Foil Spectrometer (LEFS) instruments from the Ulysses Heliosphere Instrument for 
Spectra, Composition, and Anisotropy and Low Energies (HISCALE) experiment [Hunt-Ward and Armstrong , 
2003] provide ion and electron environments at energies from a few 10’s keV to a few MeV are used to derive the 
environment. Data used in the radiation model are one hour averages of appropriate plasma moments and energetic 
particle flux interpolated to a uniform time base. 

Thermal particle populations (i=electron, proton, and helium) are modeled by Maxwellian velocity distribution 
functions 
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where the constants k b and m* are Boltzman’s constant and the mass of the particle of species i, respectively. Non- 
thermal particle populations are modeled with Kappa velocity distributions 
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Figure 1. SSRE-Rev. A Reconstructed Energy Spectra (a) The SSRE model (green) sums a Maxwellian core 
(blue) with a Kappa representation of the halo constrained by the energetic electron flux (red triangles). The 
Maxwellian halo representation is shown (red) for comparison, (b) Ion populations are modeled (green) as a sum 
of Kappa proton and helium ion distributions. Individual Maxwellian proton (dashed blue) and helium (dashed 
red) spectra are shown with the energetic proton flux (red triangles). 
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where the k>3/2. Differential flux Jj is obtained from velocity distribution functions using J i =(v 2 /m i )f 1 . 

Non-Maxwellian solar wind electron distribution functions [Montgomery et al. , 1968; Ogilvie et al ., 1971; 
Feldman et al, 1975] are often represented as a superposition of three populations. A high electron density (N c ), 
low temperature (T e ) Maxwellian core population dominates at energies <60 eV and a Maxwellian halo population 
with reduced N e and higher T e extends to energies of ~10 2 eV. An additional “super-halo” power law population 
extending from ~1 keV to ~10 2 keV exists at the highest energies [e.g., Feldman et al. , 1975; Pilipp et al , 1987; 
Larson et al. , 1996; Un et al. , 1997]. Alternatively, Kappa distribution functions can be used to represent all three 
environments [Scudder, 1992; Maksimovic et al., 1997; Meyer-Vernet and lssautier, 1998]. Figure la shows the 
combined approach used in construction of the SSRE electron differential flux spectra. A Maxwellian core 
distribution function (blue line) is added to a Kappa velocity distribution function (green line) representation for the 
halo and “super-halo” components of the electron population. Distribution functions are computed from SWOOPS 
core N c _,c T e ., c and halo N e -,h, T e . th parameters and the radial component V r of the solar wind convection velocity V c . 
Electron k parameters are obtained by adjusting the value of k until the reconstructed differential flux matches the 
LEMS and LEFS differential electron flux measurements (red triangles in Figure la). Electron populations are 
treated as isotropic at energies greater than ~10 2 keV. Asymetric electron distributions due to the field aligned 
electron strahl component and temperature anisotropies with respect to the magnetic field [Fitzenreiter et al., 1998] 
are not treated in this version of the engineering model. 

Approximatley 95% to 97% of the solar wind ions are protons with helium contributing only 3% to 5% to the ion 
number density. Other heavy ions present in the solar wind (including carbon, nitrogen, oxygen, neon, magnesium, 
silicon, and iron) only account for -1% (or less) of the total composition of solar wind ions [Ogilvie and Coplan , 
1995] and have therefore been neglected in the SSRE model environment. Proton and minor solar wind ion species 
typically exhibit non-thermal tails attached to a cold Maxwellian core ion population [Neugebauer and Snyder ; 
1966; Gloeckler et al, 1992; Chotoo et al, 1998]. SSRE differential ion flux spectra (Figure lb) are constructed 
from Kappa distribution functions for both ion species using SWOOPS proton density N H+ , proton temperature T H+ , 
and helium density N H e++. SWOOPS moments provide Tu^ and T sma n estimates of the proton temperature 
computed from two different algorithms [NSSDC/2004]. The Ti^ value, which is used for the SSRE model, is 
derived from integrating proton distribution functions in three-dimensional velocity space over all energy channels 
and angle bins. T lar ^ may yield overestimates of the proton temperature for cold solar wind flows so the choice of 
using Tiar^ provides a conservative, high flux estimate in those special conditions. SWOOPS moments do not 
include helium temperatures so the assumption T He ++ ~ T H+ is used to derive helium differential flux. Components 
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Figure 2, SSRE-Rev A 1.0 AU Flux Environments . Differential ion (top panels) and election (uottoui panda) 
flux environments (a) in the direction of the solar wind flow and (b) opposite the solar wind flow. The 
reconstructed Ulysses environments are derived from merged solar wind plasma and energetic particle observations 
scaled from the spacecraft location to 1 AU. 


of a single solar wind convection velocity V c = V r r+ V t t + V n n is reported for all ion and electron species in the 
heliocentric radial, tangential, normal (RTN) coordinate system. Only the radial component V r is used in the model 
because the tangential and normal components are typically small compared to the radial component of the solar 
wind flow. The proton k parameter is obtained by adjusting k to constrain the computed differential proton flux to 
match the measured LEMS30 and LEM120 energetic proton flux measurements (red triangles in Figure lb). The k 
parameter derived from the proton moments and energetic proton flux is simply adopted for the helium environment 
as well. 

The radius of the Ulysses spacecraft orbit varies over a range of distances from approximately 1.0 AU to 5.4 AU 
from the Sun (approximately 1.3 AU to 5.4 after Jupiter fly-by in 1992) so the solar wind samples obtained by 
Ulysses are always from radial distances from the Sun greater than the proposed Solar Polar Imager 0.5 AU or LI, 
sub-Ll missions at approximately 1.0 AU. Cold plasma data has therefore been scaled from the spacecraft location 
to the 0.5 AU Solar Polar Imager orbit or -1.0 AU LI location using l/r" scaling laws where n = 2.00 for electron 
and ion densities, n = 0.57 for ion temperatures, and n = 0.30 for electron temperatures, and solar wind velocity is 
nearly independent of radial distance beyond -0.3 AU to 0.5 AU [c.f., Hundhausen , 1972; Burlaga , 1995], 
Energetic ion, electron differential flux measurements from LEMS, LEFS instruments have not been rescaled. 
SSRE-Rev. A is applicable for both the Solar Polar Imager and the Solar Sentinels radiation environments since the 
0.5 AU environment defines the most stressing (high flux) case for both missions and 1.0 AU flux environments are 
reduced by only a factor of approximately four. 

Using the process described above, individual differential flux spectra for an appropriate radial distance from the 
Sun (0.5 AU or 1.0 AU) are computed over a range of energies 0.01 keV < E < 10 MeV in 73 energy steps for each 
set of the one hour averaged electron and ion records yielding a set of 110,572 hourly records over a 12.6 year 
period. The reconstructed Ulysses solar wind differential electron and ion flux environment shown in Figure 2 is 
the basis of the SSRE model. The format of the plot is Universal time (in decimal years) on the horizontal axis, 
particle kinetic energy on the vertical axis, and differential flux intensity indicated by the color scale where flux 
values less than the l.OxlO 2 #/cm 2 -s-sr-keV minimum threshold are set to white. 

The model includes both flux environments impacting the illuminated side of the sail computed from the radial 
component of the solar wind flow (Figure 2a) and flux environments in the direction opposite the solar wind flow 
impacting the dark side of the sail (Figure 2b). The ion environment flowing toward the Sun (impacting the dark 
side of the sail) is much reduced from that observed in the direction of the solar wind flow. However, the electron 
environment is nearly the same in both directions, particularly at energies greater than approximately 0.1 keV where 
the halo electron population in the SSRE-RevA model is nearly isotropic. 

It is important to note that the SSRE-RevA model provides only the free field charged particle environment and 
differential flux environments near the sail may be modified by the presence of photoelectrons and the Solar Sail 
potential due to spacecraft charging effects. Significant fluxes of photoelectrons emitted from the sail should only 
extend to energies of 60 eV to 100 eV and perturbations of the solar wind plasma by anticipated solar sail potentials 
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Figure 3. SSRE-Rev. A Differential Statistical Flux Environments, (a) Electron and (b) ion differential 
flux from energies of 0.1 keV to 1 MeV. Color coding indicates the 50% (black), 90% (blue), 95% (green), 

99% (orange), and maximum (red) flux environments. Differential flux to a surface in the direction 
(opposite) of the solar wind flow is indicated by the solid (dotted) lines. 

of tens of volts should only impact solar wind particles at energies on the order of tens of electron volts or less. 
Therefore, use of the environment for radiation testing at energies on the order of 1 keV and greater should be 
acceptable without considering charging effects and the photoelectron environment. The model provides fluence 
and statistical differential flux and directional flux to a surface for both the radial outward solar wind flow and flux 
moving radially inwards towards the Sun. The radial outward solar wind flow is used for characterizing the low 
energy radiation environment impact the illuminated surface of the solar sail while the radial inward solar wind flow 
provides the charged particle impact on the dark side of the solar sail. 

3.0 Discussion 

One application of the model is computation of a statistical flux spectrum from the reconstructed spectra by 
sorting differential flux within individual energy bins into monotonically increasing flux values and extracting the 
flux value which occurs N% of the way through the file where N%=50%, 90%, 95%, 99%, and 100%. The final 
100% value represents the maximum flux occurring in the record. Computing the percentile flux values for each 
energy bin yields a statistical flux as a function of energy for each of the percentile values. Figure 3 provides 
statistical electron and ion differential flux in both the direction of the solar wind (solid lines) and the direction 
opposite the solar wind (dotted line). There is little difference in the electron environment for the two directions 
since the electron thermal energy dominates the energy associated with the convective motion of the plasma. 
Differences between the two environments are only apparent at the energies less than 1 keV. The “shoulder” in the 
electron spectra from 0.1 keV to 0.5 keV is due to contributions from the Maxwellian core electron population 
during periods of high solar wind velocity and the second “shoulder” from 1-10 keV is due to the peak in the Kappa 
distribution representing the halo electron population. 

Ion environments are fundamentally anisotropic since the mean radial component of the solar wind convective 
motion is approximately 468 km/sec [Feldman et al ., 1977] which dominates thermal velocities which are typically 
only 30-60 km/sec. Multiple peaks occur in the statistical ion environment due to the effects of high and low speed 
solar wind environments and contributions from the multiple ion species. The peak with the highest flux in the 50% 
and 90% environments in particular is due to low speed solar wind protons and the secondary peak is likely a 
combination of both high speed protons and low speed helium ions. The third peak at energies of 10 keV to 20 keV 
is due to helium ions in the high speed solar wind. 

The statistical flux spectra shown in Figure 3 is a technique for evaluating mean and extreme environments 
appropriate for time dependent phenomenon associated with dose rate effects such as radiation induced conductivity 
associated with internal charging analyses. Radiation environments for estimating total ionizing dose as a function 
of depth in the solar sail material require fluence (time integrated flux) and are best obtained by integrating the 
reconstructed flux environments directly over an appropriate period of time. While it may be tempting to derive a 
mean fluence environment for use in radiation dose analyses by adopting, for example, the 50% statistical flux and 
multiplying by the appropriate length of time, this approach will underestimate the real fluence, particularly at high 
energies. 
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4.0 Summary 


A charged particle environment has been developed for use in evaluating radiation and charging issues for thin 
solar sail materials in interplanetary space. The model is based on Ulysses spacecraft observations of solar wind 
environments over a 12.6 year period. Electron and ion environments are provided for radially outward flowing 
solar wind which impacts the illuminated surface of a solar sail as well as the component of the environment 
impacting the dark side of the sail. 
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